High accuracy correction of blackbody radiation shift in an optical lattice clock 
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We have determined the frequency shift that blackbody radiation is inducing on the 5s 2 1 So - 5s5p 
3 Po clock transition in strontium. Previously its uncertainty limited the uncertainty of strontium 
lattice clocks to 1 x 10 -16 . Now the uncertainty associated to the blackbody radiation shift correction 
translates to 5 x 10~ 18 relative frequency uncertainty at room temperature. Our evaluation is based 
on a measurement of the differential dc-polarizability of the two clock states and on a modeling of 
the dynamic contribution using this value and experimental data for other atomic properties. 

PACS numbers: 32.10.Dk, 06.20.F-, 32.60.+i, 44.40.+a 



The performance of optical clocks promises a large va- 
riety of benefits. Ultimately the definition of the SI unit 
for time and frequency could be reworded to exchange 
the microwave transition as reference by an optical one. 
At present the realization of the second would be im- 
proved in accuracy by more than an order of magnitude 
[l| ; its stability could be increased even by several orders 
of magnitude 0-0] ■ Already without a changed defini- 
tion, optical clocks serve due to their accuracy to probe 
the universe in real time for temporal variations of the 
fundamental constants like the fine structure con- 

stant a or coupling of fundamental constants to forces 
like e.g. gravity [7]. In combination with their outstand- 
ing stability they are also considered to be the key instru- 
ments that could enable new measurement opportunities 
like relativistic geodesy [1, 0] with high temporal resolu- 
tion as supplement to established gravimetric techniques 
by comparing gravitational potentials directly. 

In the quest for best clocks, lattice clocks offer highest 
stability [2|-|4| but suffer from blackbody radiation shifts 
floT - ll"^ ]. Among these, strontium clocks are pursued in 
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several experiments worldwide 

systematic uncertainty now reaches the 1 x 
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FIG. 1: ac-polarizabilty of the two clock states 1 So and 3 Po as 
a function of the frequency and room temperature blackbody 
spectrum Jbbr- The resonance structure in the 3 Po polariz- 
ability is due to 115 THz 5s4d 3 Di- 3 Po transition. The magic 
wavelength lattice is at the crossing of both ac-polarizability 
curves at 369 THz. 



strated a significant reduction of uncertainty of the BBR 
shift correction by measuring the frequency sensitivity of 
the clock transition to a dc electric field Similar 
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and is thus significantly below that of primary frequency 
standards used for frequency measurements. The cor- 
rection for the shift by the blackbody radiation (BBR) 
typically dominates the uncertainty budget of the stron- 
tium clocks. Besides uncertainty due to the temperature 
measurement, at present the uncertainty of the atomic 
response to BBR at room temperature contributes with 
7 x 10~ 17 to the uncertainty. The obvious strategy to 
improve clock accuracy by reducing the environmental 
temperature requires involved apparatuses. Instead in 
most clocks accurate corrections of the BBR shift are ap- 
plied. This correction is so far only known from atomic 
structure calculations 1^, IH and lacks experimental ver- 
ification. 

As the frequency of room temperature BBR is very low, 
measurements of the difference of the dc polarizabilitics 
of the two clock states greatly improves the knowledge 
about the BBR shift. In case of an optical frequency stan- 
dard with neutral ytterbium, recent experiments demon- 



approaches have been pursued for Cs clocks [171, Il8|. 

However, in both strontium and ytterbium, more than 
the differential static polarizability Aa needs to be con- 
sidered at the level of 10~ 16 relative frequency uncer- 
tainty due to low- frequency lines (Fig. [TJ . In the BBR 
shift a dynamical correction Av^ yn (T) has to be applied 



Al/BBR = -^^ E2 )T + Al/ dyn(T), 



(1) 



where (E 2 ) T w (8.319 V/cm) 2 x (T/300 K) 4 is the mean 
squared electric field of the BBR radiation of tempera- 
ture T. To first order, Ai/dynC? 1 ) scales with T 6 . In case 
of strontium it contributes to the relative frequency shift 

-16 
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due to BBR at a level of 3 x 10" 

In this work we reduce the uncertainty of the 
BBR frequency shift coefficient of the 429 THz 
(5s5p) 3 Po - (5s 2 ) 1 So clock transition of Sr in two steps: 
First, through a measurement of the differential dc- 
polarizability Aa to provide a high precision experimen- 
tal value for the largest fraction of the shift. Second, we 



2 




Source 



correction uncertainty 



FIG. 2: Precision capacitor. Part (a) shows a sketch to scale. 
The yellow areas depict the gold plated electrodes. Part (b) 
shows the deviation of the electric field from an infinite ca- 
pacitor with U/d, when voltages ±{7/2 are applied. Part (c) 
shows the z component of the stray field Eo leaking into the 
capacitor on a plane 1 mm above the center when both elec- 
trodes are on the potential {7/2. Both fields from finite ele- 
ment calculations are normalized to U/d. 



use this observable along with several others to model 
the dynamic BBR shift, i.e. A^d yn . From the improved 
model of the response to BBR, the so far leading contri- 
bution to the uncertainty of Sr clocks can now be reduced 
by an order of magnitude to 5 x 10~ 18 at room temper- 
ature. 

To measure Aq we equipped our strontium lattice 
clock apparatus with a precision capacitor and a moving 
optical lattice setup 19] to transfer atoms from the load- 
ing to the interrogation region. Samples of few 10 4 88 Sr 
atoms at around 2 fiK are created by Zeeman slowing and 
laser cooling in a two-stage magneto-optical trap (MOT) 
20] and trapped in the horizontal one dimensional lat- 
tice with a waist radius of 65 /im and trap frequencies of 
70 kHz axial and 200 Hz radial at a trap depth of 9 /iK. 

To maintain the trap depth throughout the transport 
of atoms over several Rayleigh ranges of the focused lat- 
tice beams, the waist of the lattice beams is moved to- 
gether with the interference pattern. This is achieved 
by moving the beam shaping optics and retro-reflection 
unit with translation stages. This setup allows to move 
the atoms from the MOT region into the capacitor and, 
after interrogation by the clock laser, back out for de- 
tection by laser induced fluorescence on the (SsSp) 1 ?! - 
(5s 2 ) 1 So transition. We detect also the atom number in 
the 3 Po state by repumping them to the 1 So state to re- 
duce atom number fluctuation related noise. To avoid a 
first order Doppler shift during the interrogation of the 
clock transition due to residual motion of the stages we 
stabilize the path length between clock laser and atom 
position represented by the retro-reflection mirror of the 
lattice 
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TABLE I: Uncertainty budget for the capacitor plate separa- 
tion. 



by two optically contacted Zerodur gauge blocks of 5 mm 
height (see Fig. [5]) that was determined interferometri- 
cally using three wavelengths [22j]. The plates are par- 
tially coated with semitransparent gold layers (20 nm) 
and underlying contact layers of aluminum (5 nm) to al- 
low for an interferometric determination of the electrode 
separation before installation. The measured surface flat- 
ness of the plates leads to a correction of —153(28) nm 
(Tab. H]) at the position of the atoms. Here, we observe 
locally an angle between the electrodes of 5 /Ltrad, which 
corresponds to the measured parallelism of the gauge 
blocks' surfaces and the flatness of the electrode plates. 
After assembling the capacitor the field plate separation 
was confirmed by an independent direct interferometric 
measurement. As the plate separation was measured be- 
fore the capacitor was brought into the vacuum cham- 
ber, we also include contributions as relaxation under 
vacuum, bending of the plates under the influence of 
gravity or thermal expansion. The conductive surface 
of 28 mm x 67 mm ensures field homogeneity of better 
than 10 -6 over a large area (Fig. [5]). 

Voltage is supplied to the electrodes by two electric 
wires each, connected to the electrodes at points with 
maximum separation to also allow for measurements of 
residual resistance of the connections. Independent volt- 
ages of up to 700 V can be applied to the electrodes 
from two precision voltage supplies (Fluke 415B). The 
voltage difference is monitored by a precision voltmeter 
(HP 3458A calibrated via a Fluke 5720A calibrator to 
a Josephson voltage standard), with an averaging time 
of 1000 power line cycles (50 Hz). The electrodes are 
connected to the voltage supplies with reed relays, which 
allow for inverting or discharging the capacitor without 
interrupting the voltage measurement. 

We measure the field induced shift of the atomic reso- 
nance in 88 Sr in three field configurations (applied field, 
inverted field and no field) by interleaving three stabiliza- 
tions 



11. 23] of the clock laser to the atoms. Information 
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on the shift induced by the capacitor E cav and small bias 
The capacitor consists of two Zerodur plates, separated fields Eq is contained in the three offset frequencies 
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(i =t, I, 0) to the reference cavity of the clock laser. 



As the observed shift is proportional to ( E, 



v cap 



E, 



the dc-Stark shift due to E cap alone is given by 

fdc = (ft + - V 0, 



(2) 



provided changes in Eq are slow compared to the cycle 
times. The data of a typical experimental run are shown 
in Fig. |3] with residuals from a parabolic fit. 

The component of a small bias field Eq z along the di- 
rection of the applied field E cap is given by 



V t — V \r 00 -^cap ' E — E cap E 



(3) 



The bias field Eq z can be explained as the sum of patch 
fields independent of the applied field and stray fields. 
Patch fields can occur e.g. due to a difference in the 
work function of the gold coatings, which depends on the 
crystalline structure and can have a magnitude of up to 
500 meV [24j, |25[ . They explain the part of E 0z that is 
independent of the applied voltage (Fig. |3]). However, 
differences in the work function do not explain the volt- 
age dependence of Eq z . As this behavior may indicate 
residual differences between the voltage of the capacitor 
and the measured voltage, it was crucial to understand 
it: The system can be described by three electrodes (two 
capacitor electrodes at potentials U\ and U2 and the sur- 
rounding vacuum system on ground potential) . Any field 
in this model can be written as the sum of the antisym- 
metric configuration with electrodes at ±([/i — Uz)/2 and 
a stray field with both electrodes at the average potential 
(U1 + U2) /2. This stray field oc U1 + U2 is contributing to 
Eq z if the atoms are not probed exactly on the symmetry 
axis of the capacitor (Fig. [2] part (c)). Simulations con- 
firm that an offset of 1 mm above the capacitor's sym- 
metry plane with a shift of f mm to one side explains 
the observed dependence on the applied voltage, which 
is compatible to our effort on positioning. 

This stray field explains the voltage dependence of E 0z 
and thus Eq. [5] can be applied to determine v& c (neglect- 
ing an insignificant quadratic contribution from -Eo)- No 
significant time dependence of Eq z is observed. Along 
with measured discharging curves of the electrodes and 
measurements of the resistance of the capacitor plates 
against ground potential and between feed lines to one 
plate we conclude that the measured voltage U is deter- 
mining E cap = U/d. 

Several parabolas as in Fig. [3] (a) were measured and 
Aa was derived for each voltage. We confirmed that 
neither alternating the order of the three interrogations 
nor adding a second zero field cycle did cause any sys- 
tematic effect. No significant day-to-day variation of Aa 
was observed. To properly account for different uncer- 
tainty contributions of the voltmeter we first averaged the 
Aa values for each applied voltage. In a second step all 
these values were averaged with weights according to the 
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FIG. 3: Measured differential dc Stark shift, (a) shows the 
measured data and a fitted parabola along with the residuals, 
(b) shows the measured bias field component in ^-direction. 



uncertainty. We obtain the differential dc-polarizability 
Aa = 4.07873(11) x 10~ 39 Cm 2 /V. An influence of the 
hyperpolarizability scaling with E was investigated by 
fitting a parabola of type a + bE 2 to the Aa obtained 
in the first averaging step. The difference between a and 
the weighted average of Aa was taken as uncertainty con- 
tribution for the hyperpolarizability. A complete list of 
uncertainty contributions to Aa is given in Tab. [TXJ 

In comparison the atomic structure calculations of 
Porsev and coauthors 15] give Aa = 4.305(59) x 
10~ 39 Cm 2 /V; a difference of more than 3<r to our value. 
This is still within the typical range of deviation between 
calculated and experimental values, e.g. lifetimes in [l6j . 

As discussed in the introduction, for improving stron- 
tium lattice clocks to the 10~ 18 uncertainty range, the 
dynamic shift contribution Ai/d yn (Eq. [1]) needs to be de- 
termined. We calculated its value by first determining 
the ac-polarizabilities eei(ui) of the clock states (i = e, g) 
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fractional contribution 



(1(T 6 ) 


Separation 


15.8 


Voltage U 




Realization 


18.5 


Charging of capacitor 


5.0 


Inverting, drift 


2.0 


Residual voltage divider 


0.3 


Fringe and stray fields 


2 


Time varying bias fields 


0.3 


Hyperpolarizability 


12 


Statistical uncertainty u, U 


4.0 


Total uncertainty 


28 



TABLE II: Uncertainty budget of measurement of the differ- 
ential dc polarizability Aa. 



(Fig. [IJ using Einstein coefficients A^: 

2Jk + 1 Aki 



2J i + lu4(u4 



(4) 



and then integrating the differential ac-Stark shift over 
the Planck distribution [26[. Here u>ik denotes the angu- 
lar frequency of the transition k — > i and Ji and Jk the 
angular momenta. 

We used our measured static polarizability in combina- 
tion with other observables as magic wavelength, atomic 
lifetimes, Thomas- Reiche-Kuhn sums 27], and ac Stark 
shifts of the clock laser field [28| to adjust the Aki with 
the largest contribution in a least square fit (see Sup- 
plement). In this way we include the best experimental 
knowledge to improve A^d yn - 

In order to obtain an uncertainty estimation for Av dyn , 
we performed 2000 fit runs to synthetic datasets, i.e., sets 
of observables varied at random within their respective 
uncertainties [29(. From the Monte-Carlo simulation (see 
Supplement) we find that the dynamic BBR shifts of the 
two clock levels at To = 300 K are —150.4 mHz and 
—2.8 mHz for exited and ground state. 

In summary, the clock transition resonance frequency 
is shifted by BBR of temperature T by 



Av(T) = A^ stat — 



A^dyn 



T 



with Ai/ 8tat = 
-147.6(23) mHz. 



(5) 

-2.13023(6) Hz and Av dyn = 
Compared to the corrections used so 
far [lfj of 2.354(32) Hz at 300 K, applying our correction 
shifts the clock frequency by a fraction of —1.8 x 10 -16 . 
While this is on the order of the uncertainties of the best 
Sr lattice clocks, it is below the uncertainties of the mea- 
sured frequencies. In the course of this work a numerical 
error in the calculation of the dynamic corrections in ref- 
erence [IH was identified [jijj. This error amounts to 
—53 mHz, while this new study changes the correction 
due to the static part by 118 mHz and the dynamic part 
by 11 mHz. 



These corrections immediately apply to 88 Sr. To ob- 
tain the corrections for 87 Sr we have repeated the calcula- 
tion according to Eq. 2] using published data on isotopic 
shift and hyperfine structure 31 -3(| (see Supplement) 



and using a mass scaling of the transition dipole matrix 
elements [37|. All published transition isotope shifts be- 
tween 88 Sr and 87 Sr lie in a range of up to 150 MHz. 
Thus for the calculation for lines with unknown isotope 
shift a conservative estimate of 150 MHz was assumed. 
With these data we obtained a fractional change of Aa 
of -1.2 x 10~ 6 and of Av dyn of -3.4 x 10~ 6 . Thus the 
correction (Eq. [5]) also apply to 87 Sr well within their 
uncertainties. 

In conclusion, a measurement of the differential dc- 
polarizability of the strontium clock transition in com- 
bination with a model for the frequency response of the 
atoms to BBR enabled us to derive an improved correc- 
tion of the BBR shift and estimate of its uncertainty. 
With this study, corrected frequencies of a number of 
high precision measurements of the clock transition may 
now be calculated 0, El, El El- The corrected frequen- 
cies still agree well as the applied correction is smaller 
than typical uncertainties of realizations of the SI second 
by Cs clocks. Now the shift correction coefficient leads to 
an uncertainty of 5 x 10~ 18 at room temperature, where 
a probably achievable temperature uncertainty of 80 mK 
would add the same uncertainty. 

To further lower the uncertainty of strontium lattice 
clocks an interrogation of the atoms in a colder envi- 
ronment is a viable solution 26]. E.g. for a 1 x 10~ 18 
relative frequency uncertainty, at liquid nitrogen temper- 
ature only an uncertainty of the temperature of about 
1 K is needed, while the shift coefficient introduces a 
completely negligible uncertainty. 
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SUPPLEMENTARY MATERIAL 

This supplement covers how we obtain the dynamic 
correction and an estimate of additional corrections that 
are needed to apply our measurement on the isotope 88 Sr 
also to 87 Sr. 



FITTING PROCEDURE 

The dynamic correction Afd yn of the blackbody radia- 
tion (BBR) shift depends on the transition probabilities 
connecting the clock states to other states. In strontium 
it is mostly determined by the strength of the transi- 
tion with the longest wavelength hsAd 3 D— > 5s5p 3 Po at 
2.6 /im. Improved values of this transition amplitude as 
well as of other relevant transitions are obtained from a 
least square fit. 

In the fit we use a dataset of experimental values like 
atomic life times, the magic wavelength or the clock 
laser light shift in addition to the differential polarizabil- 
ity measured in this work. No correlation between the 
data was assumed as it originates from various unrelated 
studies. Based on Einstein A coefficients and transition 
frequencies we calculate these observables. Here we use 
the electric dipole approximation. Contributions of E2 
and Ml transitions contribute with 2.4 x 10~ 5 Hz and 
2.4 x 10 -8 Hz respectively to the room temperature black- 
body radiation shift |15| and can be neglected on the level 
of uncertainty considered here. In the fit procedure, the 
values of selected A coefficients (Type F in Tab. lIII|) were 
adjusted to achieve the best overall agreement with the 
dataset. 

For the transitions that most contribute to the polariz- 
ability, experimental line strength data is available which 
was also included in the dataset of observables (Type D 
in Tab. IIIII) . This results in a soft constraint of these free 
parameters. For all lines (Tables HVl and |V|) we have used 
experimental uncertainties to the A coefficients as given 
in the references. To the aforementioned 2.6 /jm line we 
have assigned an increased uncertainty due to the incon- 



Histogram of Monte-Carlo 
simulation (2000 runs) 
Gauss fit: -147.6(13) mHz 



sz 0.10- 




Av„ (mHz) 

dyn v ' 

FIG. 4: Probability distribution of the dynamic BBR shift 
Ai/dyn(300 K) together with a Gaussian fit to determine its 
value and the associated uncertainty. 



sistent literature values [lfj, |38[ to avoid a too strong 
constraint in the dynamic correction, which is strongly 
correlated to this value. 

The contributions of transitions to high lying ns and 
nd states and of the continuum were estimated using hy- 
drogen wave functions as in _39] . To these contributions 
an average transition energy corresponding to the energy 
difference to the ionization limit was assigned. 

Choosing other lines as free parameters in the fit 
(marked by an asterisk in Tables IIVI and [V]) did neither 
reduce the residuals significantly nor change the final re- 
sult for the dynamic correction. 

The Thomas-Reiche-Kuhn (TRK) sums for the 1 S and 
3 P systems were mainly used to avoid extreme excursions 
of single oscillator strengths that are not well constrained 
by the experimental data. Deviations from the TRK- 
rule are expected from relativistic corrections |40| . In 
the strontium spectrum relativistic corrections are on a 
fractional level of 5% of the term energies [4l| , which was 
also assumed as uncertainty in the TRK sum. 

From the magic wavelengths at 813 nm and 390 nm 
only the first one was included in the fit. The 390-nm 
wavelength depends on details in the spectral region not 
of interest here. The uncertainty of the 813-nm wave- 
length was increased to 1.5 GHz to allow for isotope shift 
of its value and corrections due to two photon transitions 
(hyperpolarizability) [42[ ■ For the static polarizability of 
the ground state we use an experimental value (43j . For 
comparison we include in our calculation an additional 
contribution of the atomic core as in (16J. This contri- 
bution cancels for all other observables, as it is the same 
for both states. 

The fitted parameters are shown in Tab. IIIII Overall 
a good agreement with the literature values is achieved, 
indicating a consistent description. 
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To estimate an uncertainty of the value of A^d yn that 
can be calculated from the fitted A coefficients, a Monte- 
Carlo simulation is performed. Artificial datasets were 
generated by randomly choosing the experimental ob- 
servables within their uncertainties. For each dataset a 
fit is performed and A^dyn is calculated. The probability 
distribution of A^d yn is shown in Fig. [4] together with a 
Gaussian fit to determine the final value and its uncer- 
tainty. Additionally we varied A coefficients not included 
in the fit by up to 15% to account for their uncertainties. 
We observed shifts of A^d yn of up to 1.0 mHz, which we 
added to the uncertainty from the Monte- Carlo simula- 
tion, resulting in Afd yn = 147.6(23) mHz. 

ISOTOPE EFFECTS 

In order to apply our results obtained with 88 Sr to op- 
tical clocks operating with the fermionic isotope 8 Sr we 
repeat our calculation, taking into account the hypcrfinc 
structure due to the nuclear spin 1 = 9/2 and the isotope 
shift. 

We rewrite the dynamic polarizability as 



47T x -> 

ati(w) = — 



u ik d; 



k 



(6) 



with transition dipole elements df k = ^^f-^fjqrr Aki/uf k 
for a transition with spontaneous emission rate A k i from 
state k to i with angular momentum J k and Jj. 

In a nonrelativistic approximation, the finite nuclear 
mass M is taken into account by replacing the electron 
mass m e in the Schrodinger equation with the reduced 
mass 



M 



(7) 



Together with the electron mass also the atomic units 
scale, e.g. for a finite nuclear mass the atomic unit of 
lengths ao is replaced by 



JL 

m c 



-ao- 



(8) 



Thus, in this approximation the dipole matrix elements 
dki scale proportionally to the reduced mass as well [37 1 . 
Comparing 88 Sr to 87 Sr, this amount amounts to 



87, 



= 1 + 7 x 10" 



(9) 



If this scaling would also apply to the transition ener- 
gies (i.e. normal mass shift oc /i), the polarizability 
would scale to the third power of the ratio of the re- 
duced masses and the dynamic correction to the fifth 
power 26]. However, in addition there are contributions 
to the isotope shift from electron correlations (specific 
mass shift) and nuclear size (volume shift). Thus, we use 
available experimental data for the transition shifts to 
calculate the difference between the two isotopes, similar 
as in high accuracy calculations of the helium [55l457| or 
lithium [58| polarizability. For lines where no data was 
available we have assumed an isotope shift of 150 MHz, 
which is larger than any observed shift. 

Where available, we have included hyperfine splittings. 
As the hyperfine interaction leaves the center of gravity 
unchanged in first order it has no influence on the dy- 
namic polarizability (Eq. [6]). Even with larg e hyperfine 
splitting on the prominent line 3 Po - 3 Di [59(, it leads to 
a correction of less than 10 -6 . Thus no hyperfine split- 
ting was used for lines where no data was available. 

From our calculation we obtain a change of the BBR 
shift of 



Ai/stat - 87 Ai/ 8tat = 2.98 /xHz, 



(10) 



or 1.3 x 10 6 of the total BBR shift. The corresponding 
change of the dynamic correction amounts to 

88 Ai/ dyn - 87 Ai/ dyn = 0.5 ^Hz, (11) 

or 3.4 x 10~ 6 of the dynamic correction. Thus on the 
level of accuracy that we have achieved, to 87 Sr and 88 Sr 
the same BBR correction can be applied. 
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Type 


property 


fit result 


literature value 


reference 


DF 


A(5s5p 1 P 1 -» 5s 2 J S ) at A = 461 nm 


-0.017 a 


1.9001(14) xlO 8 


[44] 


F 


A(4s5d 1 P 1 -)■ 5s 2 l So) at A = 231 nm 


6.6 xlO 6 s" 1 


1.7 x 10 7 s" 1 


[45] 


DF 


A(5s4d 3 Di -> 5s5p 3 P ) at 2.6 


2.67 x 10 s s _1 (-0.3 a) 


2.9(8) x 10 5 s" 1 


see text 


DF 


A(5s6s 3 Si ->• 5s5p 3 P ) at 679 nm 


7.55 xlO 6 s -1 (-2.7 a) 


8.9(5) xlO 6 s" 1 


[46J 


F 


A(5s5d 3 Di ->■ 5s5p 3 P ) at 483 nm 


3.34 x 10 7 s _1 


3.3(2) xlO 7 s" 1 


[47,48J 


F 


^4 e ff( 3 Di continuum— 5- 3 Po) at A w 316 nm 


9.71 x 10 7 s" 1 


5.6 xlO 7 s" 1 


[49] 


D 


5s 2 1 So Thomas- Reiche-Kuhn sum 


1.97 


2.0(1) 


[2JJ 


D 


5s5p 3 Po Thomas-Reiche-Kuhn sum 


1.99 


2.0(1) 


[2JJ 


D 


Aq(5s 2 % - 5s5p 3 P ) 


-9 x 10~ 5 a 


4.07873(11) x 10~ 39 Cm 2 /V 


this work 


D 


static polarizability 5s 2 1 So 


0.6 cr 


3.1(3) x 10" 39 Cm 2 /V 


[43] 


D 


dajdv at 813 nm magic wavelength 


-2.6 a 


2.2(2) x 10" 5 Hz/(£ r • MHz) 


[33J 


D 


magic wavelength near 813 nm 


3 x 10 -4 a 


368554.7(15) GHz 


[42], see text 


D 


dynamic polarizability at clock transition 


-5.4 a 


-13(2) Hz/(Wcm" 2 ) 


[28] 


u 


magic wavelength near 390 nm 


767344 GHz 


768917(18) GHz 


[50] 



TABLE III: Input data and fit results in the determination of the dynamic BBR shift. Type D denotes experimental data used 
in the least square fit, F denotes atomic parameters that are varied during the fit and u denotes experimental data that is not 
included in the fit. 



upper level 


wavelength 


A ki (s _1 ) 


reference 


5s5p 3 Pi 


689 nm 


4.69(11) x 10 4 


[51] 


5s5p 


461 nm 


1.9001(14) x 10 s 


[44J F 


5s6p 


293 nm 


1.9 x 10 6 


[45J 


5s7p x Pi 


257 nm 


5.3 x 10 6 


[45] 


5s8p x Pi 


243 nm 


1.8 x 10 7 


[45] 


5s9p 


236 nm 


1.15 x 10 7 


[45J 


5sl0p x Pi 


231 nm 


6.6 x 10 6 


[45J F 


4d5p x Pi 


243 nm 


1.7 x 10 7 


[45] 


Rydberg and cont. X V\ 


w 218 nm 


w 10 7 


estimated from [52J * 



TABLE IV: Lines connecting to the 5s 2 1 So level. Lines used in the fit are marked by F, asterisks denote lines that have been 
used as alternative free parameters. 
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upper level 


wavelength 


A ki (s -1 ) 


reference 


5s6s 3 Si 


679 nm 


8.9(5) x 10 6 


[46] 


F 


5s7s 3 Si 


433 nm 


3.1(2) x 10 6 


[47, 481 




5s8s 3 Si 


378 nm 


9.1 x 10 5 


[53] 




5s9s 3 Si 


355 nm 


5.0 x 10 5 


[531 




5sl0s 3 Si 


344 nm 


3.1 x 10 5 


[531 




5p 2 3 Pi 


474 nm 


3.9(3) x 10 7 


[47. 481 


* 


4d 2 3 Pi 


330 nm 


5.9(9) x 10 7 


[48. 541 


* 


5s4d 3 Di 


2603 nm 


2.9(8) x 10 5 


see text 


F 


5s5d 3 Di 


483 nm 


3.3(2) x 10 7 


[47, 481 


F 


5s6d 3 Di 


394 nm 


1.48 x 10 7 


[16] 




5s7d 3 Di 


363 nm 


7.9 x 10 6 


[53J 




5s8d 3 Di 


348 nm 


4.7 x 10 6 


[53] 




5s9d 3 Di 


339 nm 


3.1 x 10 6 


[53] 




Rydberg and cont. 3 Si 


m 316 nm 


1.1 x 10 7 


[39] 


* 


Rydberg and cont. 3 Di 


~ 316 nm 


5.6 x 10 7 


[39] 


F 



TABLE V: Lines connecting to the 5s5p 3 Po level. Lines used in the fit are marked by F, asterisks denote lines that have been 
used as alternative free parameters. 



